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We analyze the design of a potential replacement technology for the commercial ferrite circulators that 
are ubiquitous in contemporary quantum superconducting microwave experiments. The lossless, lumped 
element design is capable of being integrated on chip with other superconducting microwave devices, thus 
circumventing the many performance-limiting aspects of ferrite circulators. The design is based on the dy¬ 
namic modulation of DC superconducting microwave quantum interference devices (SQUIDs) that func¬ 
tion as nearly linear, tunable inductors. The connection to familiar ferrite-based circulators is a simple 
frame boost in the internal dynamics’ equation of motion. In addition to the general, schematic analysis, 
we also give an overview of many considerations necessary to achieve a practical design with a tunable cen¬ 
ter frequency in the 4-8 GHz frequency band, a bandwidth of 240 MHz, reflections at the -20 dB level, and a 
maximum signal power of approximately order 100 microwave photons per inverse bandwidth. 


I. INTRODUCTION 

With the advent of quantum information processing with 
superconducting circuits CD, the ability to route microwave 
signals without loss or added noise has become critically 
important. In particular, the operation of a digital super¬ 
conducting quantum computer will require numerous ana¬ 
log functions, such as signal amplification, feedback, and 
transduction. Moreover, it is likely that a future quan¬ 
tum information processor will employ coherent exchange 
of microwave fields among various modular components. 
Preserving quantum information in these propagating mi¬ 
crowave fields demands nearly losses and noiseless compo¬ 
nents. 

Although most passive microwave components can be 
readily fabricated with low-loss superconducting metals 
and integrated with other superconducting circuits, this is 
not true for non-reciprocal components such as isolators, 
circulators, and gyrators. Existing technology uses ferri- 
magnetic materials in intense (~0.1 T) magnetic fields to 
create the non-reciprocal behavior required to ensure the 
one-way flow of information within the network [2 -4]. Such 
magnetic devices would be quite difficult to integrate with 
superconducting circuits, which are disrupted by magnetic 
fields of 0.1 mT or less. Currently, the non-reciprocal el¬ 
ements in quantum information processing networks are 
commercially available devices connected to the rest of the 
network using meter length coaxial cables. Even if this cum¬ 
bersome arrangement were tolerable, the loss associated 
with the transition from planar circuits to coaxial cables and 
that associated with the ferrimagnetic elements themselves 
is unacceptably large 0§|8']. 

Instead of achieving non-reciprocity through the use of 
magnetic materials, one can instead use time-dependent 
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reactive elements. Long known as a general method for 
creating non-reciprocal devices fT0lfl2l . time varying reac¬ 
tances have not had much technological impact because 
ferrite elements O or active transistor [13] devices provide 
a less complex source of non-reciprocity for conventional 
electronics. But with the emergence of superconducting cir¬ 
cuit based quantum information processing, the idea of cre¬ 
ating non-reciprocal response through time varying reac¬ 
tances has returned to prominence fl4lfT7l because of the 
problems with ferrite devices and the practical absence of 
transistor technologies with quantum-limited noise perfor¬ 
mance 03 . 

Conveniently, a standard superconducting circuit ele¬ 
ment, the Josephson junction, can be operated as a time- 
variable reactance. Although the inductance of a Josephson 
junction is intrinsically non-linear, the inductance experi¬ 
enced by a small electrical signal is effectively linear, loss¬ 
less, and may be varied by also applying a larger “pump” 
current through the junction H8l. A new type of supercon¬ 
ducting non-reciprocal device exploits this effect, in which 
several oscillatory pump tones modulate the inductance 
and therefore the frequency response of several resonant 
circuits in a cyclic manner [14] [Tgl |l7l. Other novel ap¬ 
proaches to nonreciprocal circuits for superconducting net¬ 
works that don’t rely on time-variable reactances are also ac¬ 
tive areas of research [191 20]. 

In any scheme based on time-variable reactances, it is 
important that the pump tones themselves do not inter¬ 
fere with the operation of other devices in the network [17 1 . 
Furthermore the non-reciprocal device should not mod¬ 
ulate the incoming signal to create sidebands that leave 
the device. And ideally, large signals should be processed 
linearly, despite the components’ fundamental nonlinear¬ 
ities. To these ends, we introduce another concept for a 
non-reciprocal element, a four-port circulator that operates 
by only modulating the coupling rate of four itinerant mi¬ 
crowave modes to two resonant circuits in a cyclic manner. 
This four-port circulator can also be wired as a two-port gy- 
rator ED, and although not explored here, can thus also be 
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wired as a three-port circulator |2][20][22]. The symmetry of 
the circuit ensures that sidebands generated by the dynamic 
coupling are completely “erased” as they leave the device. In 
addition, the modulating pump tones do not co-propagate 
with the signal tones and oscillate at a frequency at least a 
factor of ten less than the signal frequency, and can there¬ 
fore be easily filtered out of the signal path between devices 
operating at the signal frequency. Finally, the Josephson 
junction elements are arranged into series arrays of junc¬ 
tions 1 23], an arrangement that retains the variable induc¬ 
tance but dilutes the non-linearity to approximate a linear 
time-variable inductor. 

In this article, we first analyze the general equations 
of motion of a four port circulator created through time- 
dependent coupling between the ports and the internal res¬ 
onant modes. We show that this concept is closely analo¬ 
gous to four-port ferrite circulators (where our time-varying 
reactances create a synthetic magnetic field 023) and we use 
a waveguide ferrite circulator 0 as a pedagogical touch¬ 
stone. We then introduce and analyze a circuit that re¬ 
alizes the desired circulator equations of motion. This 
non-reciprocal circuit uses a modular design comprising 
four identical subcircuits. These subcircuits are themselves 
composed of four tunable inductors arranged in a Wheat¬ 
stone bridge configuration [lT| |24] , where the tunable in¬ 
ductors are created from flux tunable SQUID arrays. In the 
final section, we analyze SQUID arrays and how they form 
flux tunable inductors. 


II. PHENOMENOLOGICAL DESCRIPTION 

Our circulator approach is conceptually related to four- 
port ferrite turnstile circulators EISj. Originally developed 
in the 1950s, these microwave circulators consist of a junc¬ 
tion of four rectangular waveguides coupled to a cylindrical 
resonator [25], Fig.[l^. As is typical of most practical circu¬ 


lators, turnstile circulators are rotationally symmetric and 
preserve the carrier frequencies of signals EQ. Nonreciproc¬ 
ity comes from a ferrite Faraday rotator in the resonator 
that rotates the polarization of signals in the cylinder. For a 
properly chosen cylinder length and rotation rate, the junc¬ 
tion acts as a four port circulator with a finite bandwidth. 
That is, when matched loads are placed on the four wave 
guide ports, microwave signals within the circulator band¬ 
width are scattered by the junction according to the scatter¬ 
ing matrix 
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where bi^t) and bi )OUt {t ) are, respectively, the complex en¬ 
velopes of a traveling wave incident on and scattered by the 
junction at port i. 

The circulator’s operation is quite modular. If there were 
no coupling between the resonator and the waveguides, the 
scattering between the waveguide ports would be reciprocal 
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with a relatively large bandwidth. This scattering matrix 
simply corresponds to the signal transfer between an input 
source driving three loads in parallel, each with the same 
impedance as the source [25]. As described in more detail 
in section [HI] when the resonator is coupled to the waveg¬ 
uides, the equation of motion for the input, output, and the 
two resonant polarization mode envelopes (assuming a sig¬ 
nal frequency at the resonator center frequency) is approxi¬ 
mately 191 (2711291 
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where a x>y are the complex envelopes of the resonant polar¬ 
ization modes (with units of (photon number) 1 / 2 ), k is the 
total energy decay rate of a resonator mode into the termi¬ 
nated waveguides, and ft is the rate of polarization mixing 
due to the Faraday rotator. Thus, the output signals, bi >out (t) 


(with units of (photon number/sec) 1/2 ), are a linear com¬ 
bination of “prompt” scattering from the input signals (as 
determined by the lower right matrix block in Eq. f3}) and 
the resonant modes (lower left matrix block). The resonant 
modes are governed by inhomogeneous, first order ordinary 
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FIG. 1: a) Schematic of a four-port ferrite turnstile circulator [3][4]|. Arrows represent the electric field polarization of the four waveguide 
and two resonator modes. The electric fields of ports 1 & 3 couple directly to resonator mode x (blue arrows), and the electric fields of 
ports 2 & 4 couple directly to resonator mode y (red arrows), b) Depiction of the circulator’s steady state operation: a traveling wave 
signal incident through waveguide port 1 induces a steady state response in the resonator modes that is rotated by angle ip relative to 
the x-mode. When ip = tt/ 4, all of the incident power is emitted out port 2. c) Concept of an alternate realization, based on Eq. (5). The 
ferrite Faraday rotator is now removed and reciprocity is now broken by the resonant cylinder mechanically rotating at rate D relative to 
the waveguide junction. 


differential equations, determined by the modes’ “internal 
dynamics” (upper left matrix block), and driving from the 
input signals (upper right matrix block). Note that in the 
absence of the rotator, H = 0, the x- and y -polarized res¬ 
onator modes couple only to waveguides 1 & 3, and 2 & 4, 
respectively. Thus, orthogonal waveguides (e.g. waveguides 
1 and 2) couple only through the prompt scattering. With 
the rotator present, excitations are “rotated” between the 
polarization modes, coupling all four waveguides through 
the resonator as well as through prompt scattering. 

Eq.® can be solved for the steady state response of the 
resonator modes to input signals. Doing so, one finds that 
the polarization of the steady state excitation in the res¬ 
onator is rotated by an angle ip = atan(2D/7c), Fig.[lJ>. When 
ip = ti/ 4, this steady state excitation couples back to each 
waveguide with equal magnitude. However, the prompt 


scattering and resonator-mediated paths carrying signals 
from {1 —> 2,2 —> 3,3 —> 4,4 —> 1} interfere constructively, 
while all other signal path interferences are completely de¬ 
structive. Thus, one finds that driving the junction on the 
cavity’s resonance, in steady state, and for ip = ti/4, the scat¬ 
tering between inputs and outputs becomes exactly Eq. (l). 

However, manipulations of Eq. f3) suggest alternative 
ways of realizing the same input-output dynamics. For ex¬ 
ample, we can rewrite these dynamics in terms of a rotating 
basis of the resonator modes 
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in which case Eq. (3) becomes 
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Inspection of Eq. (5) suggests an alternate (albeit impracti¬ 
cal) realization of a turnstile circulator. The Faraday rota¬ 
tor dynamics are absent (i.e. no explicit polarization mixing 
in the upper left hand matrix block), but the coupling be¬ 


tween the waveguides and particular resonator polarization 
modes “rotates” in time. Thus, one might imagine the same 
turnstile junction, but with the ferrite rod removed and the 
cylinder physically spinning at a rate H, Fig.fTJ:. 
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FIG. 2: a) Electrical schematic of a Wheatstone bridge-based LC resonator that serves as the basic module of our four-port circulator de¬ 
sign. Analogous to the turnstile-type structure depicted in c, the resonant mode couples to the two transmission lines 1 & 3 (characteristic 
impedance r) with a magnitude and sign that depends on 9. The resonant mode center frequency is 9 -independent, which is also true 
of c. b) Two such modules placed in parallel, with “orthogonal” bridge unbalancing, d) Four such modules, configured to emulate the 
mechanically-spinning turnstile depicted in Fig.[I]:. As mentioned in the text, the inductances in figure d are modulated linearly in section 
[II] but the inverse inductances are modulated linearly in section]!!!] 


The general structure of Eq. (5) deserves comment. In¬ 
deed, it has long been recognized that parametrically- 
modulating components in an electrical network can make 
a reciprocal network nonreciprocal fToHT2l rT4HT6l . These 
modulated components create frequency-sidebands on sig¬ 
nal carriers that serve as extra degrees of freedom to encode 
a physical location onto the signal (e.g. a signal’s “port of en¬ 
try”). Nonreciprocal networks then also require either reso¬ 
nant modes or some other means of creating delay between 
the modulated components. However, networks that don’t 
demodulate (i.e. coherently “erase”) these sidebands when 
the signal exits the network suffer either from more compli¬ 
cated downstream signal processing, or effective loss if the 
information in the sidebands is ignored. The network rep¬ 
resented in Eq. (5) achieves perfect demodulation by lim¬ 
iting the signal modulation to the coupling between reso¬ 
nant modes (with no internal loss) and the itinerant fields. 
Thus, when signals pass through the resonant modes they 
are modulated exactly twice, as they enter and as they exit 
the resonant modes. In general, only 0 th and 2 nd order har¬ 
monics of the modulation frequency can appear in the out¬ 
put signal, and it is easy to design these modulations such 
the 2 nd order harmonics completely cancel out. In contrast, 
many other, related proposals modulate the coupling be¬ 
tween resonant modes or between itinerant fields, or some 
combination of all three modulation types 1TQUT21 [T4UT61 . 
Unfortunately, these types of modulation create sidebands 
at all harmonic orders, which are more difficult to cancel or 
to filter out. 

Our approach emulates the dynamics in Eq. (5), but as 
a lumped element, superconducting microwave network 
with simulated “spinning” of two resonant modes relative 
to four transmission line ports. Our network can be built up 
piecewise, starting with the network depicted in Fig. [2^. In 
the limit of small inductors (relative to the port impedance) 
the two transmission lines 1 and 3 are effectively shorted to¬ 
gether through the bridge of inductors over most frequen¬ 
cies. However, the bridge also presents a total inductance 
of l to the capacitor, forming a resonator with center fre¬ 
quency l/VTc. For cos [6) / 0, the bridge unbalances, cou¬ 


pling the resonant mode to the transmission lines with a 
magnitude and sign that depends on 6 |l'7i,|24]. The anal¬ 
ogous turnstile-type structure is depicted in Fig. |2]: where 
a two-port rectangular waveguide couples to a rectangu¬ 
lar box resonator ||25]. Assuming that only one polarization 
mode of the box is near resonant with signals applied at the 
ports, the magnitude and phase of the coupling depend on 
the orientation angle 6 of the box relative to the waveguide. 
Adding a second bridge resonator in parallel with the first, 
but now with the inductors imbalanced as ±sin(9 ), Fig. If. 
would be analogous to replacing the rectangular box res¬ 
onator in Fig. [2]: with a rotationally symmetric, cylindrical 
resonator (like that in Fig. [!]:). Finally, the full circulator de¬ 
picted in Fig. [TJ: is emulated by adding two more transmis¬ 
sion line ports 2 & 4 that couple to the same two resonant 
modes, but do so “orthogonally” relative to ports 1 & 3 via 
two more bridges. If the unbalancing of the bridges can be 
varied sinusoidally in time, i.e. make 6 = ftt as in Fig. If, 
then the two resonant modes couple and uncouple from the 
four transmission lines in a coordinated fashion that sim¬ 
ulates the mechanical spinning of the cylinder in the four 
port circulator depicted in Fig.fTJ: and in Eq. (5). The trick, of 
course, is to find an experimentally convenient way to both 
realize dynamically tunable inductors in a superconducting 
microwave circuit and tune them naturally in such a highly 
coordinated fashion. 

Our tunable inductors are dc superconducting quantum 
interference devices (SQUIDs), which are two Josephson 
junctions connected in parallel l26l . For currents much 
smaller than the SQUID critical current, the SQUID acts as 
an inductor with inductance 


ls = 
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( 6 ) 


(assuming identical junctions and negligible geometric in¬ 
ductance) where I s is the SQUID critical current, J 0 the 
junction critical current, </? 0 the reduced flux quantum, and 
$ the magnetic flux threading the SQUID loop. Replac¬ 
ing the four inductors in Fig. [2^ with four SQUIDs real¬ 
izes a bridge of dynamically tunable inductors. Further- 
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FIG. 3: Schematic of a dc SQUID-based realization of the tunable 
bridge network. In the right hand figure, the wire configuration 
(black lines) is critical as it determines the magnetic flux through 
the five depicted loops (four SQUID loops, one loop of the bridge 
itself. The four small circles represent unimportant network ter¬ 
minals.). The total magnetic flux $ through each SQUID loop is 
the sum of the flux from a uniform background field ($ s ) and from 
a magnetic flux control current that flows vertically through the 
“twisted” bridge ($ A ). Ideally, the net magnetic flux through the 
bridge loop is always zero. 

more, twisting the bridge into a figure-eight layout gives 
us a convenient way to achieve the required unbalancing 
of the bridge inductances. We can set the inductance of 
each SQUID appropriately by applying a constant and uni¬ 
form magnetic field and a smaller, time-dependent, and 
gradiometic magnetic field, as in Fig. [5] Note, too, that this 
configuration only induces the desired screening currents 
within each SQUID loop. Undesired screening currents 
between SQUIDs are not induced (canceling unwanted 
screening currents in a network like Fig. [2 ] i is also impor¬ 
tant, but outside the scope of this article). Thus, in princi¬ 
ple, all of the inductive, dynamically unbalanced bridges in 
Fig.[2]l, could be realized with 16 SQUIDs, a uniform back¬ 
ground magnetic field and two control wires carrying cur¬ 
rent oscillating as cos(Qt) and sin(fit). In section 0 we 
will also address considerations such as increasing the net¬ 
work’s saturation power by replacing single SQUIDs with ar¬ 
rays and the secant- rather than linear-dependence of in¬ 
ductance on magnetic flux. 

III. CIRCUIT ANALYSIS 

We now analyze the circuit depicted in Fig. [2]l. Analysis 
will use both the frequency domain, lumped element ap¬ 
proach common in microwave engineering [25), and a time- 


domain approach, which shows how Eqs. [5} approximate 
the circuit’s dynamics. 

The essential sub-network in this circulator is the two 
port [25j] Wheatstone bridge-type (24) network depicted in 
Fig.[4^. In contrast to the networks depicted in Fig. [2] in this 
section we will exclusively consider inductances that vary 
as /(I ± d(t, 0)) _1 (/ some inductance, and d is a real func¬ 
tion of time t, and an angle 6), as in Fig. [4^. This is a more 
natural, simplified model for the dependence of SQUID in¬ 
ductances on magnetic flux, Eq. 0. Straightforward circuit 
analysis [25] gives us the constitutive equations for this net¬ 
work 
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oo Vi(z)dz is the time integral of the volt¬ 
age across port i and is called the “branch flux” across this 
port l30l . and It is the current entering at port i. Through¬ 
out this article, we will use the convention that inductance 
is defined as the ratio of the branch flux across and current 
through a network branch H8ll26l (which for time-varying 
or nonlinear inductances is different from the more famil¬ 
iar definition of inductance as the ratio of voltage and the 
time derivative of current). From Eq. (7), we learn that 
the unbalancing of the bridge in Fig. [5^ gives an output 
reluctance (the current response at one port per flux ap¬ 
plied at the other; i.e. inverse inductance) proportional to 
d(t,0), but the input reluctance (the current response at 
one port per flux applied at the same port) is independent 
of d(t,Q). Such clean separation between input and output 
reluctances is highly attractive when d{t, 6) is proportional 
to an experimentally-convenient tunable parameter, such 
as the control current in Fig. [3] Our network is designed to 
exploit this separation to realize time-variable coupling be¬ 
tween resonant modes and itinerant fields without modu¬ 
lating the intra-resonator dynamics or prompt field scatter¬ 
ing. This bridge is the repeated module that helps us realize 
Eq. (5). 

Next, Fig.|4jD takes four copies of this bridge, places them 
in a ring, and defines six ports {1,2,3,4 ,q,p}. Note that 
ports 1-4 are each defined between a network node and 
ground, while ports q and p are each defined by two net¬ 
work nodes. If we let d(t,6) = ecos(fit + 9) (0 < e < 1 
), straightforward circuit analysis then relates port branch 
fluxes and currents: 
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While the structure of Eq. (8) is evocative of Eq. they 


are not equivalent. For example, Eq. 0 contains no reso- 
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FIG. 4: a) The inductance bridge as a network with two ports, 1 & q. b) Four inductance bridges configured as a six port network. Attaching 
a capacitor across ports q & p and transmission lines across ports 1-4 realizes the circulator network depicted in Fig.[2]i. 


nant dynamics yet and its dynamics are reciprocal H0ll25l . 
Nonetheless, we already have the critical structure that only 
the couplings between ports 1-4 and ports q & p are vari¬ 
able in time. Analysis of Eq. (8) simplifies by going into a 
“left-right/even-odd” basis for the port 1-4 variables and a 
rotating, circular basis for the q and p port variables: 
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and similarly for the current variables {j being the imag¬ 
inary unit, adopting the electrical engineering conven¬ 
tion l38l ). In these bases, the dynamics of Eq. (8) separate 
into: 


they turn the reciprocal, four-port “odd” network Eq. GD 
into a nonreciprocal, two port network. 


A. Frequency-domain analysis 

Taking into account the capacitors and writing the odd 
network dynamics now in the frequency domain, we find: 
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The time-dependent, reciprocal network of Fig. [5] d be¬ 
comes nonreciprocal when we make modes q and p reso¬ 
nant by placing a capacitor (capacitance c) across each of 
these ports. Doing so, these ports gain a fixed relation be¬ 
tween current and branch flux 

d 2 

~ C dt 2 = (13) 

or equivalently in the circular basis 


TjnJ 0±(O = (14) 

The capacitors have no effect on the “even” dynamics in 
Eq. (lT), but in the following subsection we will show that 



where square brackets [•] indicate a frequency domain vari¬ 
able. From Eq. GD, one finds that a current applied at ei¬ 
ther the “left” or “right” odd ports (i.e. taking 7/ )0 / 0 or 
I r>0 / 0) induces a resonant response in the 0± variables 
when co = \J (4- e 2 )/2/c ± D = co 0 ± D, respectively. For 
D = 0, the orthogonal coupling between the odd modes and 
the resonant modes means that the left and right odd modes 
are uncoupled over all drive frequencies. But, for ft ± 0 and 
odd driving at the frequency co 0 , the 0+ and 0_ modes re¬ 
spond as equal and opposite capacitive and inductive reac¬ 
tances in parallel, which open a transmission window be¬ 
tween the left and right odd ports (2]. For instance, elimi¬ 
nating 0± in Eq. (15), setting co = co 0 , and expanding to first 
order in D we find that 
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In this configuration, signal transfer through the network is 
nonreciprocal: voltage applied at the right odd port will in¬ 
duce a positive current at the left odd port, but voltage ap¬ 
plied at the left odd port will induce a negative current at the 
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left odd port (note that jcoo times a branch flux is a voltage). 
A network with the constituent relations given in Eq. fl6) is 
known as a “gyrator,” the most basic building block of non¬ 
reciprocity in electrical network theory, with “gyration resis¬ 
tance” e 2 /16cf7 l[2T1 . 

Further analysis is aided by attaching transmission lines 
to ports 1-4. With the capacitors across ports p and q 
and transmission lines across ports 1-4, Fig. becomes 
equivalent to Fig. [2jl. The voltages and currents at ports 
1-4 are now decomposed into traveling wave voltages and 
currents incident upon (“in”) and scattered by (“out”) the 
network through the transmission lines [25] (27; [29]. Our 
goal is now to relate these input and output waves. While 
such input-output analysis is convenient theoretically when 
studying four-port circulator networks [39], measuring scat¬ 
tering parameters is also more experimentally convenient 
in microwave networks than direct measurements of port 
voltages and currents. 

If the only constraints on the input and output fields are 
those imposed by our network (e.g. the transmission lines 
are of effectively infinite length) then I9l l25ll27l 

00) = 0,in0 ) T 0,out0 ) 

d 

“TT0*0) = M]in(f)T Vi out(f) 

dt 

V li, inO) — M]in0) 

— Y 0, outO) = 14,out 0) (17) 

where the traveling wave variables above are evaluated 
at the port positions and r is the characteristic (real) 
impedance of each transmission line. Standard network 
analysis [25] then tells us that if the port branch fluxes and 
currents are related by 

j coY i j [co] 0 j [co] = h [co] ( 18 ) 

where ¥ is known as the admittance matrix, it then follows 
that 


HoutM = [(1 + r¥[co]) 

= S tJ [<o]V ]M [{o], (19) 

where 1 is the identity and provided that the matrix inverse 
exists. S is often referred to as a scattering matrix. 

For example, using Eq. m one finds that when 17 = 17 () = 
e 2 /16cr 


signal picks up a n phase shift traveling from the right odd 
to left odd port and no phase shift traveling from left odd 
to right odd. Similarly, using the even network relations in 
Eq-dT}. one finds that to zeroth order in co 0 l/r 


Vl,e, outlet)] 


'or 

Vl, e , in[^>o] 

Tr,e,outlet)] 


1 0 

kr,e,in[^o] 


( 21 ) 


Combining Eqs. (20} and (21} and putting the input and out¬ 
put voltage signals back into the port 1-4 basis, we find that 
this network acts as an ideal, four-port circulator, Eq. (l), for 
input voltage waves at frequency co 0 . We also mention that 
as our network is operable as a two port gyrator, it is there¬ 
fore also operable as a three-port circulator by combining 
ports 3 and 4 in Fig.[2]d into a single port (i.e. by merging the 
transmission lines 3 and 4 together) l20ll22l . 


B. Time-domain analysis 

This standard scattering matrix approach, though, ob¬ 
scures the dynamics that guided our intuition: a simulated 
“rotation” of the resonant modes relative to the itinerant 
fields (e.g. Fig. [ 2 ]:). To that end, we reanalyze the dy¬ 
namics of Eq. (Tl) and Eq. CD in the time-domain, and 
without eliminating the 0± dynamics. In quantum optics, 
analogous models are known as “input-output” (10) mod¬ 
els [9 2"7-- 29]. Tike the scattering matrix models of the pre¬ 
vious paragraph, 10 models relate incident and scattered 
electromagnetic waves and are thus natural models for res¬ 
onant microwave networks. But unlike scattering matrix 
models, they are time-domain, first order ordinary differ¬ 
ential equations. 10 models approximate resonant circuit 
dynamics, but may also be directly compared with, for ex¬ 
ample, optical systems that don’t have lumped element rep¬ 
resentations. On the other hand, the circuit analysis in the 
previous section does not require approximations beyond 
the lump element assumptions, and can be useful in con¬ 
sidering how 10 approximations break down. 

In this case, our 10 model assumes that the network vari¬ 
ables have solutions of the form | 9 ]27 28 

0*0) = Wi{t)e jC0dt + c.c. 

V/,inO) = CL>oV~crv i)in {t)e jc0dt +c.c. 

KoutO) = oJoV~crVi )OUt {t)e^ OJdt +c.c. 

cod = &>o + A. (22) 


V\ } o, outt^o] 


' 0 -1 ' 

Vl,0,m[o>o\ 



1 0 

Vr,o,i n[^o] 


In other words, when the gyrator resistance equals r, the 
odd network is matched: signals are fully transmitted by the 
network, but the transmission is nonreciprocal. A voltage 


for a drive frequency detuned by A from the center fre¬ 
quency co 0 . Then, combining Eqs. (TT), (15}, and (TT), mak¬ 
ing a slowly varying envelope approximation and solving to 
lowest order in (o 0 l/r , |A|/&> 0 , and \Q \/&> 0 (i.e. assuming 
that the rate of variation of </?*, Vi t m, and Vi >out is much less 
than cc>d) gives 
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h- 


-W 


Vl,o,out{t ) 
Vr,o, out(T ) 


(j( A-n)+|) 0 

, v'T -/I 

0 -fj(A + f2)+|) 


-N T ~h T 

-1 0 

-vT V2 

0 -1 


<p+{t) 


Vl,e, out(T) 


0 1' 

Vl,e,m[t ) 

Vr,e, out(f) 


1 0 

Vr,e,in{t ) 


Vr,o, in(0 


where we have defined k = e 2 /8cr. We call Eq. {23} an 
10 model of the network depicted in Fig. |2ji. Putting 
this model back into the {1,2,3,4 ,q,p} basis then gives 
us the four port circulator model Eq. {5} with the iden¬ 
tifications {CLq, CLp, b (in/out) t ^2,(in/out )> ^3, (in/out) > ^4, (in/out)} = 
{(fq> (fp> ^i, (in/out) > ^ 2 , (in/out )> V 3 , (in/out )> ^4, (in/out) }• A natural in¬ 
terpretation of the 10 model Eq. {5} is that while the cen¬ 
ter frequency and total damping of each resonator mode 
are constant in time, the resonant modes couple to (and are 
damped by) each transmission line in a sinusoidally “rotat¬ 
ing” pattern. While the rotating coupling between resonant 
and port degrees of freedom is apparent in Fig.[2jl, the time- 
independent intra-resonator dynamics are less obvious. 

Finally, we calculate the performance of the circulator 
in this 10 representation in order to compare it to the fre¬ 


quency domain circuit analysis presented in section III A 
For linear 10 models of the form 


] 


' A 

[B ' 

x{t) 

_ youtO) _ 


. € 

D _ 

Tin(f) _ 


(24) 


with time-independent matrices {A, B,C, D}, a scattering 
matrix-like representation of the steady state response of 
the output amplitude vector y out to a constant input ampli¬ 
tude vector y 0 ut is 


you, = (D-CA-^)?" (25) 


in general. However, as is typical for 10 models of lossless, 
gain-less, linear systems, Eqs. {23} have a particular form 
where B = -C+D, A = —j Q - IC+C, and D is unitary and 
Q is Hermitian (+ signifying the matrix adjoint) l29l . In this 
case, in analogy with Eq. {no, it can be shown that [ 32] 

y™ t = D 1/2 (1 + Y io r 1 (1 — Y, 0 ) D 1/2 >f n s 
Y i0 = 2yD 1/2 (CQ- 1 C + ) _1 D- 1/2 (26) 

(assuming the matrix inverses exist). Using Eq. {26} and the 
first Eq. in {23} , we find an admittance-like matrix for the 10 
model of the odd modes 


-2 


J io,odd : 


/ A —ft 
n y'A 


(27) 


And thus, in analogy with Eq. {20}, we find for A = 0 and 
ft = k/2 (i.e. for cod = (o 0 and ft = ft 0 ) this 10 acts like an 
ideal gyrator 


" v ss 
l,e, out 


' 0 

-1 ' 


r vf s . 1 

l,e, m 

V ss . 


1 

0 


v ss ■ 

r,e,out 





r,e, m 


(28) 


If tuned optimally {ft = ft 0 ) and for A / 0, the FWHM band¬ 
width of t he odd an d full circulator network models are 
a/2 k and \j 2(a/3 — 1)tc, respectively. Thus circulator band¬ 
width is strongly constrained by the modulation amplitude 
of the inductors, k oc e 2 y as will be discussed more in sec¬ 
tion [iv| We also note that a/2 k = 2V2ft 0 , which is V2 times 
the frequency splitting of the internal modes for a matched 
network. A simple proportionality between bandwidth and 
the frequency splitting of internal modes is typical of many 
non-reciprocal networks [2]. 


C. Simulation 

It is worth comparing the theoretical scattering response 
of the approximate 10 model to the full lumped-element 
scattering matrix (which is exact up to the lumped element 
assumption). In the 10 model, the scattering response cor¬ 
responds to an ideal four port circulator when A = 0 and ft = 
k/2. Alternatively, one can calculate S for the full lumped 
element network depicted in Fig. [2}hbut without the ap¬ 
proximations made to produce EqsT{2Q) -{21}. Choosing pa¬ 
rameters compatible with superconducting microwave cir¬ 
cuits mm l = \ nH, c = 2 pF, r = 50H, co = coo = 2nx 6.16 
GHz, ft = ft 0 = 2n x 99 MHz, and e — 1 we find: 


ISnl 2 = IS22I 2 = |S 33 | 2 = |S 44 | 2 = 0.002 
IS21I 2 = IS32I 2 = IS43I 2 = |S 14 | 2 = 0.995 
IS31 1 2 = IS42I 2 = IS13 1 2 = IS24I 2 = 0.002 
IS41I 2 = IS12I 2 = IS23 1 2 = IS34I 2 = 0.000 (29) 

That is, 99.5% of the incident power at this frequency is 
routed correctly for these parameters and 0.2% is reflected, 
assuming ideal lumped element components and ideal in¬ 
ductor modulation (section [Tv| will focus on the many con¬ 
siderations that go into identifying a practical set of pa¬ 
rameters). The network’s response to off-resonant drives, 
A = cod ~ / 0, is depicted in Fig. [ 5 ] for the various 

models. For these parameters, the circulator bandwidth 
is \] 2(a/3 - 1)k/2tz = 241 MHz in the IO model response, 
which is depicted in Fig. [5^. We see a very similar response 
using the lumped element circuit model in Fig. [ 5 J 3 . The 
greatest difference in this response is that that it is no longer 
perfectly symmetric about A = 0, especially outside of the 
circulator bandwidth. 

The discrepancy between the lumped element and IO 
models appears to be largely due to nonidealities in the 
even network transmission Eq. {2T}: perfect circulation 
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a) 10 network 


b) Lumped element network 
with luo/r = 0.4 


c) Lumped element network 
with luj 0 /r = 0.8 





- |Sn | 2 

- |S 2 i | 2 

- |S 3 i | 2 

- |S 4 i | 2 


FIG. 5: Simulation of the frequency domain response of various models of the circulator, as described in the text. Curves of the other 
scattering parameters are identical to these, with a cyclic permutation of the port indices, confirming proper circulation. 


breaks down when the even network ceases to look like two 
ports shorted together. Thus the ideal even network re¬ 
sponse occurs in the limit of small co 0 l/r, and for these pa¬ 
rameters coq l/r = 0.39. If instead we choose the parameters 
l = 1 nH, c = 1 pF, e = 1/a/2 (co 0 = 2n x 6.7 GHz), the cor¬ 
responding 10 model still has the same co 0f k , and H 0 , but 
now coqI/t = 0.84. Evaluating the lumped element scatter¬ 
ing parameters for these values now results in a response 
that deviates further from the ideal IO model, with a 97.8% 
maximum power transmission, 1% power reflected, and a 
more asymmetric response about A = 0, Fig. [ 5 ]:. 


IV. CONSIDERATIONS FOR A SQUID-BASED REALIZATION 

To better ground the above analysis in the reality of su¬ 
perconducting microwave circuits, we now outline some 
considerations that go into designing dc SQUID networks 
that approximate linear inductors. Circulator performance 
should be compared against commercial circulators and the 
needs of contemporary and future experiments with super¬ 
conducting quantum microwave circuits. Desirable charac¬ 
teristics for our circulator include an input-output response 
that is close to Eq. (l) on resonance, a wide bandwidth, and 
high power handling. For example, there is a trade-off in 
this particular circuit design (Fig. [2]l) in that ideal input- 
output response occurs in the limit of small (Ool/r, subsec¬ 
tion [III]:, but the bandwidth of the network also decreases 
in this same limit ( k/co 0 ~ lcoo/r). This trade-off is inde¬ 
pendent of any component nonidealities and arises at the 
schematic level of this particular design. We also note that 
competition between bandwidth and the optimal perfor¬ 
mance at a single frequency appears in any practical circuit 
design. The physics of dc SQUIDs also require trade offs in 
a good design. Although modifications to this circuit will 
change the particulars, the considerations outlined below 
appear in any superconducting microwave network based 
on dynamically-modulated SQUIDs I6l l23ll26ll3ll . 

We first address the schematic-level trade-off. One must 
assume that r, the transmission line impedance, is prac¬ 
tically constrained to 50 D and that we will target designs 
in which coo/2n, the center frequency of operation, is in 


the usual 4-8 GHz band for quantum superconducting mi¬ 
crowave networks. In optimizing the ratio co^l/r, this leaves 
only l as the flexible parameter. As described in sub¬ 
section [1111 :, if network bandwidth is held constant, but 
this ratio increases, ideal circulator operation will decrease. 
For example, finite co 0 l/r produces back reflections (e.g. 
ISnl 2 / 0 in Fig. [ 5 }. While it is possible to include addi¬ 
tional, lossless matching networks at each port to reduce 
back reflections over some frequency range, the Bode-Fano 
criterion suggests that a nonvanishing co^l/r fundamen¬ 
tally limits excellent network matching over broad band- 
widths |25 33 . This further supports the intuition gained 
from Fig.[5]that without a radical redesign of the network, a 
larger / decreases the performance achievable over a given 
bandwidth of operation. 

To roughly match the performance of commercial ferrite 
circulators, we can allow ourselves nonidealities in the sig¬ 
nal routing at the -20 dB level. We see such performance 
in Fig. [5]:, in which l = 1 nH. Inductances of this scale 
are achievable in SQUID networks fabricated using optical 
lithography I35ll36l . Although the performance depicted in 
Fig. [5J3 is better than in c, the assumed parameters are not 
realistic. The baseline impedance is a reasonable l = 0.5 nH, 
but Fig. [ 5)3 assumes 6 = 1, which implies that the SQUIDs 
periodically achieve infinite inductance (the feasibility of 
e = l/V2 will be discussed below). Infinite inductance is ob¬ 
viously unphysical, but the actual limitations require some 
discussion of SQUID physics. 

A single dc SQUID is formed by two Josephson junctions 
connected in parallel by superconducting wires, forming 
a loop, Fig. [6^. For simplicity, we assume the junctions 
are identical and junction capacitance is negligible, and we 
ignore the geometric inductance of the wires for the mo¬ 
ment. Then, Faraday’s law and the Josephson relations al¬ 
low us to relate the voltage across (Vi) and current flowing 
between (Ji) two leads connected on opposite sides of a 
single SQUID loop in the presence of an external magnetic 
field [261. Assuming null initial conditions, we identify a sin¬ 
gle SQUID’s effective inductance L\ as the ratio between the 

/ £ 

-00 W T )^ T ) i n SQUID 
leads. SQUID inductance is nonlinear, as this ratio depends 
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a) 


b) I N 


1 





FIG. 6: a) Schematic of a dc SQUID with J 0 critical current foseph- 
son junctions, l g geometric inductance in each branch, and a mag¬ 
netic flux b) A series array of N such SQUIDs. 


on 7i, and we find that to third order in I\/I s 


(pi (o = 

Li(f) sa 

hit) = 


L\{t)Ii{t ), with 

wo (i+!(/(* vw*)) 2 


)■ 


= S ec f 5^1 

Is{t) 21 0 \2ifo J 


(30) 

(31) 

(32) 


where (p 0 = H/2e is the reduced flux quantum, I s is the 
SQUID critical current (the maximum current flowing be¬ 
tween these two leads that the junction can support in its 
superconducting state), / 0 is the critical current of a single 
junction, and <!>(/) is what the magnetic field flux through 
the loop would be in the absence of screening currents. 
Thus, for small currents, a SQUID acts like a linear induc¬ 
tor with a magnetic flux-dependent inductance. The rela¬ 
tive magnitude of the linear and nonlinear inductances, 1 / 6 , 
is fixed: a larger inductor requires a SQUID that saturates at 
smaller currents. 

Now consider an array of N identical SQUIDs with their 
leads connected in series, Fig. [ 6 J 3 . In this case, the current 
flowing through each, and the voltage across each are iden¬ 
tical. Thus, the relationship between the branch flux across 
and the current through the entire array is (p N = L N I N , 
where L N = NL\. Qualitatively, adding SQUIDs in series al¬ 
lows us to increase linear inductance without decreasing the 
saturation current. Quantitatively, to third order in h/I a we 
have that |23 34| 


L N (t) t* + (33) 

WO = TT7V Ia(t) = Is(t)/N, (34) 


from which we see that for a fixed linear array inductance l a , 
the nonlinear inductance scales as N~ 2 . Less nonlinearity 
means higher saturation powers in analog signal process¬ 
ing. 

The SQUIDs’ nonlinearity limits the circulator’s power 
handling. As the SQUID inductance increases with cur¬ 
rent, center frequencies of microwave resonators contain¬ 
ing SQUIDs typically decrease as the power of incident sig¬ 
nals increase. The approximation of a SQUID as a linear 
inductor tends to break down when this center frequency 
shift is comparable to the resonator bandwidth. In the case 
of high quality resonators, it is convenient to reparameter¬ 
ize the nonlinear inductance in Eqs. (33 and (33) as an ef¬ 
fective “Kerr constant” that gives this center frequency shift 
per microwave photon stored in the resonator (i.e. stored 
microwave energy in units of ficoo) [23 [34j. Up to a factor of 
order unity (depending on the resonator construction) this 
Kerr constant is 


Kf 


tKOp 

1*1 


■ocAT 


(35) 


where the sign of K indicates the direction of the center fre¬ 
quency shift. For example, using parameters compatible 
with Fig. [ 5 J 3 , co 0 = 271 x 6.16 GHz, l a = 1 nH, and assuming 
an array of IV = 20 SQUIDs compatible with junctions pro¬ 
duced by optical lithography in the NIST NbAlOxNb trilayer 
process [35, 36], I s = 6.6/iA, one finds K ^ -2 n x 580 kHz. 
With this Kerr constant, and assuming the resonator band¬ 
width depicted in Fig. [ 5 ]:, one would expect the nonlinear 
reactance of the SQUIDs to become noticeable with of or¬ 
der 100 photons stored in the resonator (because 241 MHz 
= {2ti)~ 1 \K\ x 415 photons). A comparable design based on 
single SQUIDs would require I s = 0.33/iA and would ex¬ 
hibit nonlinearities well below 10 stored photons. For quan¬ 
tum superconducting microwave experiments, the satura¬ 
tion power of commercial ferrite circulators is effectively 
infinite. While experiments involving one or two super¬ 
conducting qubits often operate at signal powers that do 
not saturate routers and amplifiers networks with single- 
SQUID-scale saturation currents, higher saturations powers 
are more critical in quantum networks that employ more 
weakly coupled systems (e.g. mechanical oscillators and 
magnetic spins), applications outside of quantum informa¬ 
tion (e.g. astrophysical detectors), and in current and future 
quantum networks with more signal multiplexing. 

Optimizing our circuit design according to the prescrip¬ 
tions thus far, one might conclude that the best strategy is 
to employ SQUID arrays with l a & 1 nH (with no external 
flux bias) as our tunable inductors, and let N —> 00 , I 0 —> 0 
to minimize nonlinear effects. This is imprudent for several 
reason. One is a larger circuit footprint, which adds stray 
geometrical inductance throughout the circuit. Other rea¬ 
sons are due to the geometric inductance of the SQUID loop 
wires themselves, l g in Fig. [6^, which we’ve ignored thus 
far. This SQUID loop inductance limits both the minimum 
^ 5 ,min and maximum / 5)max achievable through flux tuning 
of a SQUID array. Both constraints worsen as (fo/Iolg de¬ 
creases |26|. With minimum and maximum array indue- 
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tances in general one finds that 


component yields effective inductances that vary as 


kl< 


rf — 1 
r] 2 + \* 


*7 — ^s, max/ Is, n 


(36) 


In practice, r] seems to be limited to about 4 in supercon¬ 
ducting microwave circuits with of order 10-100 SQUIDs in 
an array [23], which limits |e| < 0.88. Even longer arrays 
would be even less tunable. Constraints on the maximum e 
ultimately limit circulator bandwidth, which is proportional 
to e 2 (section[ni|3). The simulation depicted in Fig. [5]:, with 
e = 1/a/ 2, is consistent with this tunability constraint, and 
achieves a bandwidth of 241 MHz. Commercial circulators 
used in contemporary superconducting microwave experi¬ 
ments have bandwidth of a few GHz, but the vast majority 
of this bandwidth is unused. Most circuits involved in rout¬ 
ing quantum microwave signals in contemporary experi¬ 
ments operate at bandwidths of order 10 MHz or smaller [3- 
[7l H7l[23l[3ll . Moreover, the center frequency of our network 
may be tuned broadly by varying the uniform background 
field, as depicted in Fig. [5] Thus, while not as broadband as 
commercial options, this network should have an achiev¬ 
able bandwidth more than sufficient for most applications 
in superconducting microwave networks. 


1. Input-output response with SQUID-like inductance 
modulation 

We saw in section [In] that proper periodic modulation of 
the inductances according to Fig.|2j4 causes the network to 
scatter signals without creating any sideband excitations. 
Using a numerical, time-domain, MATFAB Simscape sim¬ 
ulation of the circuit depicted in Fig. [2]l, network perfor¬ 
mance in this ideal case is depicted in Fig. [7^. Using the 
same parameters as in Fig. [ 5 ]:, r = 50 fi, Z = 1 nH, c = 1 pF, 
e = l/\/2, U = U 0 = 2n x 99 MHz, a continuous-wave, in¬ 
put current drive at fequency cod = 2n x 6.66 GHz is applied 
to port 1 suddenly at time t= 0. The top plot in Fig. [7^ de¬ 
picts the resulting current measured in the output signals 
for the first 5 ns of the simulation (after separating input 
from output signals in the simulation). Initially, the output 
current is equally distributed among the four ports, but af¬ 
ter a few ns nearly all of the output current is measured at 
port 2. The bottom plot gives the normalized output cur¬ 
rent power spectrum for a 250 ns simulation (i.e. well into 
steady state). The spectra are sharply peaked at cod , with the 
peak of the Ii >out > h,out> and I 4j0Ut signals at least 20 dB below 
that of J 2 , 0 ut- While the shoulders are slightly fat, there are 
no visible sidebands. Driving the inputs of the other ports 
produces responses identical to Fig. [7^, up to a cyclical per¬ 
mutation of the port indices. 

Unfortunately, it is not simple to tune a SQUID’s induc¬ 
tance exactly as 

Z'(f) = /(l±ecos(Ht)) _1 or /'(£) = Z(1 =besin(Ut)) _1 (37) 

as in Fig. ii- Because the linear inductance of a SQUID 
or array of SQUIDs goes as Eq. applying a simple back¬ 
ground magnetic flux that consists of a static and sinusoidal 


W) = 


V 0 


2/ 0 


cos 


V 2^0; 


4>(Z) = $x±$ACOs(Ht) or $(t) = ±$ A sin(Ut). (38) 

It is possible to make l s exactly equal to l' s in the limit 
(fo/2Io —> 0, $x/2(^ 0 —> tt/ 2. However, this limit corresponds 
to an infinitely flux tunable Josephson inductance, which is 
unphysical. Instead, it is more realistic to set y?o/2Io = 2/, 
where l is the desired l in Eq. {37}, and set $ s /2(/? 0 = tt/3 so 
that l s = l' s when $a = 0. One then chooses $a such l s ^l' s , 
given the desired e in Eq. {37) . 

We can now redo the simulation in Fig. [7k, but with a 
variable inductance model based on Eq. {38} rather than 
the ideal Eq. {37}. Choosing parameter values that closely 
emulate the simulation in Fig. [7^, r = 50 H, (p 0 /2I 0 = 2 nH, 
$x/2(/?o = tt/3, c = 1 pF, $a/2(/?o = 0.38, H = 2tt x 90 MHz, 
cod =2n x 6.63 GHz (parameters found through trial-and- 
error optimization), produces the analogous plots in Fig.[7]3. 
The time-domain response is nearly identical to the ideal 
case. The most noticeable difference is that the / 3)0U t re¬ 
sponse is about the same as Ji j0U t and / 4)OUt . The differences 
between the steady state power spectra are more obvious. 
The output signals now exhibit sidebands, which are caused 
by the slight differences between Eqs. {38} and {37} for these 
parameters. Happily, though, the largest of these sidebands 
are nearly -50 dB c and are the 3rd harmonics of H (i.e. are at 
00 d ± 4H). The lower harmonics of H are almost completely 
suppressed by the symmetry of the circuit, which is remark¬ 
able considering the significant, 4 >a/2^o = 0.38, modula¬ 
tions of the SQUID-like inductances. 

A quantitative analysis of the spectrum observed in 
Fig. [7}} is beyond the scope of this article, but will be con¬ 
sidered in a companion publication (which will also contain 
a quantized model of the circulator) l37l . This is done by 
first deriving an Hamiltonian describing the circulator in¬ 
cluding its input and output ports. In a second step, Yurke’s 
approach to input-output formalism is used to compute the 
scattering properties of the circuit l27l . In contrast to the 
analysis in section III B| following Yurke allows us to obtain 
exact, time-domain equations of motions for the circuit’s 
degrees of freedom. Using theses results, the spectrum ob¬ 
served in Fig.[7j3 can be understood quantitatively. 


V. CONCLUSION 

We have introduced and analyzed a novel design for 
a four-port circulator based on time-varying inductances 
(Fig. |2]i) that has the potential to replace the lossy and non- 
integrable commercial ferrite circulators ubiquitous in su¬ 
perconducting quantum microwave experiments. The ba¬ 
sic module is a dynamically-unbalanced bridge network 
of four SQUIDs (or SQUID arrays) [3], depicted in Fig. [ 3 ] 
These bridges dynamically modulate the coupling between 
four transmission lines and two resonant modes, and the 
construction ensures that this coupling is the only aspect 
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Ideal Inductance Tuning 
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b) SQUID-like Tuning 
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FIG. 7: Numerical network simulations, a) Network output response in time and frequency, driving port 1, given parameters from Fig. & 
and the ideal inductance modulation given in Fig.[2]i and Eq. {37}. b) The same network simulation, but where the linear inductances are 
modulated in a “SQUID-like" manner, i.e. using Eq. {38}, as described in the text. 


of the dynamics that is modulated. The construction also 
ensures that when the inductances are tuned perfectly, the 
network’s output signals are not complicated by frequency, 
phase, or amplitude modulation. This critical feature is 
greatly aided by limiting the dynamic modulation to the 
coupling between the resonant modes and the transmission 
lines. As a result, the resonant modes act as if they are “spin¬ 
ning” relative to the transmission lines or, equivalently, as if 
they are coupled by a synthetic magnetic field, Fig. [I] 

The design has a number of attractive features such as 
modulation pump tones detuned from the signal by more 
than a decade, a broadly tunable center frequency, a tun¬ 
able bandwidth, high saturation power, and a lumped ele¬ 
ment and modular construction. We have given an overview 
of many of the considerations needed to realize a design ca¬ 
pable of achieving a 240 MHz bandwidth and a saturation 
power of order 100 microwave photons per inverse band¬ 
width with superconducting microwave technology. The 
design of a practical electromagnetic circulator is a dif¬ 
ficult and important problem and different technological 


constraints have yielded a wide variety of solutions and 
potential solutions over the decades, of which Refs. |2 - 
IU EMzl HMD are but a small sample. Because of this 
long tradition, we analyzed this network from several differ¬ 
ent perspectives, through analogy to familiar, ferrite-based 
circulator designs, through lumped element circuit analy¬ 
sis, through approximate input-output models (appropri¬ 
ate for distributed RF or optical systems), in the time and 
frequency domain, analytically, and numerically. We ex¬ 
pect that the present design will find near-term use context 
of quantum superconducting microwave networks, but we 
also hope that this diverse analysis will help inspire an even 
greater diversity of future solutions and perspectives in even 
broader contexts. 
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